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A PENNY-SHAPED CRACK IN A 
FILAMENT-REINFORCED MATRIX* 

I. THE FILAMENT MODEL 

by 

F, Erdogan and A. H. Pacella** 
Lehigh University, Bethlehem, Pa. 


Abstract. This study deals with the elastostatic problem of a 
penny-shaped crack in an elastic matrix which is reinforced by 
filaments or fibers perpendi cu 1 ar to the plane of the crack. An 
elastic filament model is developed in the first paper. The 
second paper considers the application of the model to the penny- 
shaped crack problem in which the filaments of finite length are 
symmetrically distributed around the crack. The reinforcement 
problem for the cracked matrix with elastic fibers of different 
diameter, modulus, and relative location is considered in the 
third paper. Since the primary interest is in the application 
of the results to studies relating to the fracture of fiber or^ 
filament-reinforced composites and reinforced concrete, the main 
emphasis of the study will be on the evaluation of the stress 
intensity factor along the periphery of the crack, the stresses 
in the filaments or fibers, and the interface shear between the 
matrix and the filaments or fibers. 


1 . INTRODUCTION 


The primary objective of this series of papers is to develop 
a technique by which, with a reasonable amount of computational 
effort, one may obtain the solution of the three-dimensional 
elasticity problem for a matrix containing a penny-shaped crack 
and reinforced by elastic filaments or fibers perpendicular to 
the plane of the crack. Basically, the problem is one of 

* This work was supported by The National Science Foundation 
under the Grant GK-11977 and by the National Aeronautics and 
Space Administration under the Grant NGR-39-007-011. 

**Present Address: Bettis Atomic Power Lab., Westinghouse 

Electric Corp., West Mifflin, Pa. 
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Interaction between a macroscopic crack and filaments or fibers 
in a composite medium. The problem finds its practical applica- 
tions in the fracture studies of fiber or filament reinforced 
composites and in reinforced concrete. Even though the existence 
of such a crack has an effect on the vibration characteristics, 
the stiffness, and other mechanical properties of the material, 
its main importance lies in the reduction it causes in the frac- 
ture resistance of the structure. Hence, our primary attention 
will be concentrated on the calculation of such quantities as the 
distribution of the stress intensity factor along the periphery 
of the crack, the filament-matrix shear stress and the maximum 
tensile stress in the filaments. 

This first paper in the series will be devoted to the devel- 
opment of a model for an elastic filament imbedded into an elastic 
matrix. The main requirements expected of the model are a suffi- 
ciently accurate representation of the filament, and its applica- 
bility to the interaction problems involving a cracked elastic 
continuum with mu 1 ti -f i 1 amen t reinforcements. For a sparsely 
reinforced matrix in which the interaction between the perturbed 
stress fields of the isolated filaments and the crack is negli- 
gible, the solution given in [1] for an ellipsoidal inclusionin 
an infinite matrix may be quite satisfactory provided the filament 
ends are rounded and there is no excessive concentration of inter- 
face shear. However, since the filaments are usually cylindrical 
with sharp edges and since the technique described in [1] cannot 
readily be expanded to interaction problems, the ellipsoidal 
inclusion model of [1] is not suitable for the problem under 
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consideration. A somewhat more appropriate model for the present 
purpose would be that described in [2]. The model described in 
[2] would give sufficiently accurate results for the tensile 
stress in the filament and for the stiffening effect on the crack. 
However, its representation of the filament-matrix contact 
stresses would not be sufficiently accurate. Partly for this 
reason and partly for reasons of convenience in solving the 
resulting integral equations, in this paper a somewhat different 
model will be developed. The technique is based on a direct 
generalization of the notions discussed in [3] and [4], and will 
be described in the next section. Some numerical examples will 
then be given and the results will be compared with those obtained 
from using the methods of [1] and [2]. 

2. SOLUTION OF THE GENERAL INCLUSION PROBLEM 

Consider the three-dimensional inclusion problem shown in 

Figure 1. Let the homogeneous, isotropic elastic domains (the 

inclusions) which are bounded by nonintersecting smooth surfaces 

S|^, (k = 1 ,2, . . . ,m) be perfectly bonded to the surrounding elastic 

medium (the matrix). Let the bounding surface S^ of be 

n ^ 

subjected to surface tractions T^. (where S^ may be finite or 
infinite). Let the elastic constants of Dj^ be (k = 

0,1,. ..,m). The problem may be formulated by writing the field 
equations for separately with the boundary conditions on 
and stress and displacement continuity conditions on S-j,...,S^. 
This, however, requires the solution of an elasticity problem for 
the simple domains as well as for the multiply-connected 

domain D^. The problem may also be considered as that of a 
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Figure 1. General inclusion geometry. 




simply-connected non homogeneous domain In which the elastic 
constants have jump discontinuities along the surfaces 
This formulation requires the solution of a problem in which the 
field equations have discontinuous coefficients. Aside from the 
special case discussed in [1], neither one of these solutions is 
tractable. However, it can be shown that the problem may be 
reduced to the solution of a system of integral equations pro- 
vided the Poisson's ratios of the elastic domains , D-| , . . . , 
are assumed to be equal. For certain geometries these integral 
equations may be solved numerically without any difficulty. 


Let D and S be the union of inclusion domains and their 

boundaries, respectively, i.e., 

m m 

D = I D|. , S = I S|^ 

1 1 


( 1 ) 


Let u-|,U 2 ,U 2 be the components of the displacement vector in the 
nonhomogeneous medium (D^ + D+S), In the absence of body forces, 
the elastostatic boundary value problem may be formulated as 


+ = 0 , (Xj e (D+S + Dp)), 

°ij ^ ’ f’”’ =1*2.3), 

where the discontinuous elastic constants are given by 

y ~ Up > ^ ~ s ^^j^ ^r* r~0,l,...,m), 

0 n ^ . 

n. is the outward normal and T- is the traction vector on 

J ' 

In (2) -(4) as well as in the rest of this paper the usual 


( 2 ) 

(3) 

(4) 

(5) 

^0* 

summa- 




tion convention is used. Let 



V = + Au , A = + AX ; (6) 

Au = - Uq , AA = A^ - A^ , (Xj e D^, r=1,...,m), (7) 

where Ay and AA are assumed to be nonzero. With (6), (2) may be 
expressed as 


y u • 
^0 1 


, JJ 


+ 




J1 


+ [Ayu- •• + (AA+ Ay)u. ..] = 0 

1 » j J J > J ' 


[x. e (Djj + D + S) ; i = 1 ,2,3} . (8) 


In (8) the quantity in brackets may be considered as a body force 
vector which, due to the discontinuous nature of the coefficients 
A and y, is expected to be discontinuous across and at the bound- 
ary S. Let r\ = ri(x.)» (x,j € S) be the distance along the normal 
n measured (in outward direction) from the surface S. We may 
then define 


Ayu ■ . . + ( AA + Ay ) u . - • = ^ ( x . ) , ( x - € D + S ) 

ijjj jjj* ‘J J 




pyxj) , 


(Xj 6 D) , 


( 9 . a , b ) 


•Tyxj)6(n - 0) , 


(Xj e s). 


where, in general , is not equal to the boundary value FT. 
Using (9), (8) may now be expressed as 

(Xj « (D(, + D+ S), i= 1,2,3). 


( 10 ) 


f 

On the other hand if we let F-(x-) = F. , (x.£ D„) and 

1 J 1 J r ' 

Ti(Xj) " '■'i ’ ^ equivalent to 

^^r^i,jj ^ = 0 > ^ Dp)* 
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(n -a-c) 


r r _n 


(Xj e s^) . 


* “j,i^ ^ ^^“k,k^-j • 

( r = 1 , . . . » m ; i , j = 1 , 2 , 3 ) . 


(10) with (3) and (4), and (11) give the formulation of m+1 elas- 
ticity problems for the simply-connected homogeneous domains 
(Dq+$+D) ,D^ , . , . ,Dj^. In addition to displacement components 
u^ =u? in (Dq+D+S) and uV in the auxiliary inclusions with 

the elastic constants Ay and AA^, (r=l,...,m; i= 1,2,3), the 

r 

equations contain the unknown functions and T^. , (r= l,...,m; 
i= 1,2,3). The additional equations to account for these un- 
knowns may be obtained by considering the fact that the displace- 
ments on S are continuous and u- given by (11) are identical to 
u. given by (9), (8), or (10), namely 

Ui(x.) = u°(Xj) , (Xj € D^), 

u^{xy = u°(xy , (x. e s_), (r = l . ,m; i = 1,2.3). 

(12. a, b) 

0 r r n ^ 

Formally, 9m+3 unknown functions u^^, u,j, , T,j may be obtained 

from 9m+3 equations given by (10), (11. a), and (12) under the 

boundary conditions (3) (with (4)) and (11. b) (with (11. c)). If 

the Green's functions for the domains (D^+D+S) (i.e., the simply- 

connected region bounded by S^) and S^, (r = l,...,m) are known, 

this system of equations may easily be replaced by a system of 

r n*^ 

6m integral equations for the unknown functions F- and , 

(i = 1 ,2,3; r = 1 , . . . ,m) . 


Let us now consider the field equations for the part of 
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the homogeneous medium (D^+D+S) in which u^. (Xj) = u^. (xj). This 
may be obtained from (10) as 


1 


“i.jJ ^ ° ’ 


(x-^ Dj r l,...)m) i»j l»2j3)* 

«J ’ 

r o_ r 


(13) 


subject to the boundary conditions that “ ^ij'^j* ^ ^r^ 

where is the limit of the stress component obtained from the 
J 

solution of (10) as x. approaches the boundary from inside. 

vJ 

From (13), (11. a) and (12. b) it is easily seen that 
From (14) it then follows that 

fT = 0 . if ^ ^r’ = 1 » • • • i = 1 ,2,3), (15) 

where v is the Poisson's ratio of the elastic region D , 
s ^ 

(s = 0,1 , . . . ,m) . Thus, with the assumption that , 

(r=l,...,m) the formulation of the problem may be considerably 
simplified and may be summarized as 


V?,jj * + J, T? = 0 ■ 

(Xj € (Dg+D+S)), 

" “Li’ " VLk«iJ^"; = (X. € S„), 06. a. b) 

^^“i,jJ ^ "^^r)“Lji = ° ’ ^ °r’’ 


[Au^(u:^. . uj^.) . = - T. 


(Xj e S^). (17. a, b) 


u,.(Xj) = ii?(Xj), (Xj.eS^). (r=l,...,m; 1,j= 1,2,3). (18) 
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Again, if the Green's functions for the regions for a 

concentrated stress vector on the boundary and for (D^+D+S) for 
a concentrated internal body force are known, (18) directly gives 
a system of two-dimensional integral equations for the unknown 
functions , (i = l,2,3; r=l,...,m). 


In the study of the mechanics of composite materials an 


important quantity of interest is the magnitude of the contact 
stresses on the interfaces Once the layers of body 

nr 

forces T- , (r=l,.««,m) are determined, the contact stresses 

may easily be obtained from the equilibrium considerations along 

the boundaries $ 1 ,.,.,$^^. Let be the components of 

m 

the stress vector on the internal surface S, (5=2 S^) having 

the normal (n^.)» (i,j = 1,2,3). Let x^ be the components of the 

contact stress vector on the interface having the outward 
r 

normal n^- . It is clear that 


t;.(Xj) = T^(Xj) = °;k(Xj)n' (Xj) , (Xj e s'). 


0 + 


(19) 


where the superscripts + and - refer to the boundary values of 
the related quantities as the surface is approached from out- 
side (the positive side) and from inside (the negative side), 

respectively. The equilibrium considerations for the homogeneous 

r\^ 

region (D^+D+S) subjected to the layer of body forces T^. on S^, 

n ® 

(r=l,...,m) and surface tractions T. on require that 

T?+ - T°- + t" =0. (Xj € S^; i = 1.2.3; r= (20) 

Also, from the solution of the problem for the simply-connected 


domain (see (17)) we have 

r_ r, r _n'^ 

= °ij"j - T- • 


(xj^e s^). 
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On the other hand, If , ( r = 1 , . . - »m) , using the equality 

of the displacements in D^, (xj ^ D^, r=l,.**>ni), 

from the stress displacement relations it may easily be shown 


that 


a!v ’ (X^e D^. r= (22) 


or 


1 


^0 


AU, 




(23) 


Thus, from (19), (20), (21), and (23) the components of the con- 
tact stress vector on S may be obtained as 


;.(xj) = X^"(x.) = - (Xj.) , 


0 + 


(x-€ S; i,j=l>2,3; r l,...,m). (24) 

0 * 


Once the problem for the m+1 simpl y«connected domains is 
solved, noting that the displacement components u,|, (i = 1,2,3) 
in the actual inclusion with elastic constants and are 
given by 

Ui(xp = u°(Xj) = ur(x.) , 

(x-^ D ; i,j = l»2,3; r = l,...,m), (25) 

xJ * 

and (because of = v^, (r=l,...,m)) 



the stresses in the actual inclusion may be expressed as 




'V(“i,j + “j,i> ^ V“k,k«ij 




(x^ e D^; r = 1 , 
i , j = 1,2,3), 


. ,m; 


( 27 ) 
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where a?., (i,j== 1,2,3) are the stress components in the matrix 

. r 

(Dq+D+S) which has the elastic constants X^, and cr.j, 
(i,j=l»2,3) are the stress components in the auxiliary inclusion 
with the elastic constants Au^,, AX^, (r= l,,..,m). 


Note 1 . The results found in this section remain valid for the 
plane and axisymmetric problems, with the additional simplifica- 
tion that for vp = vq, (r=l,.«.,m) the resulting integral equa- 
tions for the unknown functions TO^, (i =1,2; r= l,..*,m) would 
be one-dimensional, ^ ^ 


Note 2 . In the corresponding "antiplane shear" problem 


(28) 


and the results found in this section regarding the vanishing of 
the body forces remain valid without any restriction on the 
elastic constants ur» ( r = 0 , 1 , , , . ,m) . In this case (16-18) and 
(24) (with (25)) give the exact solution. In this problem too 
the interface may be represented by enclosed plane curve and 
the resulting integral equations (for , r=l,*,,,m) are one- 
dimensional, the arc length measured along 5^ being the variable. 


3. THE FILAMENT MODEL 

Let the filament be represented by a cylindrical inclusion 
of length 2c, radius r^, and the elastic constants E^, v^. Let 
the elastic constants of the surrounding matrix be E, v. It is 
assumed that 

(a) = v; 

(b) the dimensions of the matrix are large in comparison 
with c; 

(c) the external load is the traction applied to 

the matrix away from and parallel to the filament; and 

(d) the length of the filament, 2c, is large in comparison 

with its diameter Thus, the following basic relations for 

the infinite medium may be used in deriving the Green's functions 
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for the matrix [5] : 
“1 


(i = 1 ,2.3), : (29) 


= !(x, - x.)^ , 
1 

. _ (1 + v){3 - 4v) 

^ " 8ttE{ 1 - v) 


B = 


1 + V 

8ttE( 1 - Vr 


(30) 


where u,j , (i = 1,2,3) are the components of the displacement 
vector at the point x. due to the concentrated body forces, 
acting at the point Tj , (j =1,2,3), and x,j and Xj refer to the 
rectangular coordinates* If we deal with an axisymmetric problem 
in which, referred to the cylindrical coordinates r,0,z, the body 
forces R,@,Z are distributed over a ring r=r^, 0<e<2'ir, z = t 
in such a way that 0 = 0 and R and Z are independent of 6, 
integrating over the ring, from (29) the displacement components 
at a point (^"'^q» O_<0<2tt, z) may be obtained as 

u.(r^,z) = K^^(z,t)R + Kip(z,t)Z , 


r ' 0 


u^(r^,z) = K^-|(z,t)R + Kpp(z,t)Z 


Ug = 0 , 


‘21 


22 


11 


(z,t) = [2r„ + illxilizli!. ] K(k) 


(31 . a-c) 


- — [p„ + 7 ^- (2r^ + (t-z)^)]E{k) 

r o p o 

0 


K^2(z.t) = - K2,(z,t) = - [K(k) - E(k)] 


<22(2. t) = 


4Ar 


[K(k) + yE(k)] , 


( 32 . a-c) 
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Y = i ^ ^ k = (33) 

where K(k) and E(k) are complete elliptic integrals of the first 

and the second kind, respectively. A list of integrals used in 

the derivation of the kernels K.*, (i,j=l,2) may be found in 

• 0 

Appendix A, Similar expressions may be obtained for 
u^(r,+c), u^(r,+c) due to R,Z at (r^,t), 

u^(r,+c), u^(r,+c) due to R,Z at (s,+c), and 

u^(r^,z), u^(r^,z) due to R,Z at (s,+c), 

where 0£(r,s)<r^, ~c<(z,t)<c. 

The filament model developed in this section will be used to 
study the stress state around the leading edge of a penny-shaped 
crack in the matrix located in the z=0 plane. Since r^ << c and 
since the body forces R are locally self-equilibrating, the direct 
effect of R on the stress intensity factors along the crack 
periphery will be negligible. However, since the integral equa- 
tions in R and Z will be coupled, the effect of R through Z may 
not be negligible. The first example discussed in this section 
will be devoted to study the effect of neglecting R on Z. For 
the sake of simplicity and in order to consider an extreme case, 
it will be assumed that the inclusion is rigid and the end 
effects are negligible. Thus, if the uniaxial tension 
is the external load applied to the matrix away from the inclu- 
sion region (see the insert in Figure 2), from (31) the integral 
equations of the problem may be expressed by writing the displace- 
ment components along (r^r^, -c<z<c) equal to zero as follows: 
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12.0f- 


Fi gure 


10.0H 


i^cO.35 

C/to=10 

CAo- 20 



2. Radial and axial contact stresses for a rigid filament. 




vr^a^ c 

u.(r-.,z) = ^ + / [K^i(z.t)R(t) + K^ 2 (z,t)Z(t)]dt = 0 , 

" ^ “C 

Uz(ro,z) = ^ + / [K2^(z,t)R(t) + K22(z , t ) Z { t ) ]dt = 0 , 

(-C < z < c) , (34) 

where K.., (i,j =1,2) are given by (32). A close examination of 

* \J 

the kernels around z = t would indicate that K^-j and K 22 have 
logarithmic singularities. This may be seen by observing that 
at z = t E(k) is finite and for small values of |t-z| we have 
the following asymptotic relation: 

K(k) = - 1og|t-z| + log 4p^ 

2 

+ 1 [- log|t-z| + log 4pg - 1] + ... (35) 

PQ 

Since the system (34) is of the first kind it is equivalent to a 
system of singular integral equations. In order to extract the 
correct behavior of the solution, it would be simpler to cast the 
system in the standard form with Cauchy-type singularities by 
formally differentiating the equations. Thus, separating the 
singular parts of the kernels, (34) becomes 

/ (z,t)R(t)dt 

-c -C 

C 

- k -j 2 ( 2 ) 1 ^ ( i' ) ^ i' = 0 , 

^0 -c 

^ / [log|t-z| + k, 2 (z.t)]R(t)dt + / 

0 -c -c 

c 4ir(l-v)a 

+ / k 22 (z*t)Z(t)dt = - (i+v)(3-4vy ’ (-c < z < c), (36) 


- 13 - 



where y = 1/(3 -4v) and the bounded functions ^ , (i , j = 1 ,2) are 


given by 
^11 


^''o , aK(k) 1 ' . ^'^0 ' ^0 + (Uy)(t-z)^ 3Kikl 
^ 3Z ■ t-z' Po(t-z) r^Pg 3z 


+ K(k)[^ - 

o Y' r> 


0^0 
2 


r p 
0^0 


. J_ i|^ [p + ^ {Zr-l + (t-z)2}] 

r^ 8z ‘*'^0 pp 0 


t-z 

•^oPo 


E(k)[l + 2y - ^ {2rl + (t-z)^}] 
PQ 


k,2(z.t) = k2^(z,t) = 2?^ - Po^l°9l*-" 


- [K(k) + log|t-z|] + -^ E(k) 

PQ 0 

+ lilzii [K(k) - E(k)] + . 


k (z t) = ^ - T^] + ' K(k) 

■^22' ’ ^ P L t-z-* p (t-z) 


2r^-P„ 2r„(t-z) 

'o 


2»‘„Y 2 c/t\ 2r v(t-z) 

+ — — E(k) , 


Pq 


(37.a-c) 


iMJli = ilii - K(k) , 

o Z "C “ Z Z 


^ [E(k) - K(k)] . 
^0 


2 _ 


4r 


4r^ + (t-z)^ 

n ' 


(38) 


Referring to [6] it may be shown that the solutions of (36), 
R and Z, have integrable singularities at +c , the index of the 
system is k= 1, and hence the general solution will contain two 
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arbitrary constants. On the other hand (36) states that the 
z-deri va ti ves of the displacements and rather than and 
are zero along (r= r^, -c < z <c). Thus, the solution of (36) 
must satisfy a set of si ngl e-val uedness conditions which may then 
be used to determine the arbitrary constants resulting from the 
general solution. These conditions may be expressed by fixing 
u^ and u^ at any point along the line of integration, say, for 
example at z = 0, giving 

u^Cr^j.O) = 0 , ° ’ (39. a, b) 

where the expressions for u^ and u^ are given by (34). 

Considering the symmetry of the problem, the solution of the 
system of singular integral equations (36) subject to the condi- 
tions (39) is of the following form [6]: 

R(z) = F(z)(c^ - , Z(z) = G(z)(c^ - (40. a, b) 

where F(z) = F(-z) and G(z) = -G(-z) are bounded functions which 
may easily be obtained numerically (e.g., [7,8]). Some numerical 
results obtained from (36) are shown in Figures 2-4. Figure 2 
gives R(z) and Z(z) for v= 0.35, (c/r^) = 10 and (c/r^) = 20. 

Figure 3 shows the effect of the Poisson’s ratio (v=0.2 and 
v=0.35 used in the Figure roughly correspond to a glass and an 
epoxy matrix, respectively). From the viewpoint of this study 
aiming to simplify the filament model the important result is 
shown in Figure 4. Here the body force Z(z) (which, in this case 
is also the contact stress) is given as obtained from (36) with 
and without neglecting R(z). It appears that for the practical 
range of c/r^ ratios the effect of neglecting R on Z will be 
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Figure 3. The effect of Poisson's ratio on the contact 
stresses for a rigid filament. 


Is ^ 




Figure 4. The effect of the radial body force on the axial 
contact stress in a rigid filament. 


IS2> 


negligible. Hence, for the remainder of this study the radial 
component R(z) of the body force will be neglected. 


In the case of the elastic filament, in order to simplify 
the resulting system of integral equations, in addition to 
neglecting the radial component R of the body force, it will be 
assumed thaf the body force Z(r,+c), (0£r< r^) at the ends is 
uniformly distributed. Again, the effect of this assumption will 
be local and will be negligible on the stresses in the matrix in 
z=0 plane (and hence, on the stress intensity factor along the 
leading edge of the crack when the crack problem is considered). 
Thus, the unknown quantities will be the distributed body force 
Z(rQ,z), (“C<z<c) and the constant Z(r,c) = p = -Z(r,-c). 

These quantities will be determined from the integral equation 
and the algebraic equation obtained by matching the displacements 
of the matrix and the auxiliary filament (with elastic constants 
t = E.-E and v) along the surface (r=r^, -c<z<c) and at an 

CL T ^ 

appropriate point at the end z=c (which will be selected as 
r=0, z=c). 


Due to the large 1 ength-to-di ameter ratio c/r^, the filament 

will be approximated by a one-dimensional body subjected to body 

forces - -J- Z(z) distributed uniformly over the cross-section 
“o 

{0<r<r^i z) and the end tractions -p distributed again uniformly 
over the ends z = +c. Thus, the displacement in the filament may 
be expressed as 


u^^(z) = -u 


fz 


(-z) 





z c 


/ dt / Z(n)dn 
0 t 


> 


(0 < z < c) , (41 ) 
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( 42 ) 


Q 

' ■ r ~ ~' f ' / tz(t)dt] . 

“"f ■ '" 0 0 

Evaluating now the displacements u^(rQ,z), (-c<z<c) and u^(0,c) 
in the matrix due to the body forces Z(r^,t), (-c<t><c) and 
Z(r,c) = -Z(r,-c) = p, (0£r<r^), and the traction at infinity 
we obta i n 

Uz(ro,z) = + / K 22 (z,t)Z(t)dt - Ap[M^(z) - N^(z) 

-c 


+ y{(c+z)^M2(z) - (c-z)^N2(z)n 


u (O.c) = + 2Trr A / (c + ^ ^ ^ ^ ^ 2-.1 /2 

^ ^ ° -c r^ + (c-tr [r^+ c-tr] ' 

0 0 

+ 2 TTAC 2 PC , (43^35b) 


where 


Y = 


3 - 4v * 


^ = (1+v)(3 - 4v) 
8ttE ( 1 -v) 


= 2 + - (4 + 


r'/c2)l/2 . 


4y( 


1 


(4 




- V 


(44) 


M.(z) = / 

0 

^0 

Ni(z) = / 

0 


rdr 


rdr 


ZlT 

/ 

0 

2tt 

/ 

0 


d0 

[p2 ^ 2rr cos0 + ( c + z ) ^ ] ^ ^ ^ ^ ^ ^ ^ 

de 

[p2 + ^2 _ 2rr cose + { c-z ) ^ ^ 

(i - 1,2,...), (45. a, b) 


and K 22 given by (32. c). Again, differentiating (41) and 
(43. a) and using (42) and (43. b), from the conditions of conti- 
nuity of the displacements we find 
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/ + / k22(z.t)Z(t)dt + ^ / Z(t)dt 

2C E 

+ f [1^2(2) + , (-c<z<c), 

c 2 C,E c 

/ k^(t)Z(t)dt + y,-7E — 1 rr ! tZ{t)dt 
-c 0 ' ‘"f '0 


C E 

+ pc(C2 + ^ Ve ) * 

t 


(46. a, b) 


where 

^ _ 47r(l-v) 

" (l +vT( 3 - 4v) * 

the kernel k 22 ( 2 ,t) is given by (37. c) and 

ki(t) = (1 + ^), 

' [r2 + (c-t)2]l/2 (3 - 4v)[r2 + (c-t)2] 

= n - + (c-z)N2(z)] 

+ 3[(c+z)^M3(z) + (c-z)^N3(z)] . (47. a, b) 

The functions and N.| appearing in (47. b) are defined by (45). 
Some simplifications for the evaluation of these functions may 
be found in Appendix B. 


The integral equation (46. a) and the algebraic equation 
(46. b) determine the function Z(z) and the constant p. Noting 
that at z = 0 K 22 is an even function of t, from (43. a), (45) 

and (41) it is seen that the si ngl e-val uedness condition 

^ will be automatically satisfied provided 
the solution of (46) is restricted to a class of odd functions 
(as required by the symmetry of the problem), i.e., Z(t) == -Z(-t), 
(-c<t<c), and Z(r,c) = p = -Z(r,-c), (0£r<rQ). The numerical 
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solution of (46) may again be obtained in a s trai ghtforward 
way [7,8]. 

Once Z(z) and p are obtained all the desired field quanti- 
ties may be evaluated in terms of definite integrals having the 
related Green's functions as kernels and 1 and p as density 

functions. In fracture studies, of particular interest are the 
+ 

contact shear along the filament-matrix interface and 

the axial stress filament. The general expres- 

sion for the contact stress is given by (24), which in this case 
becomes 

The general expression for the stresses in the filament is given 
by ( 27 ) , namely 

Of2z(i".z) = |z| c), (49) 

where is the stress in the matrix due to the external load? 
Z(z), and p, and a,,, is the axial stress in the auxiliary 

0 a Z Z 

filament which has the elastic constants E = Ep - E and v. cr,, 

CL f 2.Z 

appearing in (49) may be obtained by adding to the stress 
component evaluated from (29) and the related stress- 

displacement relations. Here, since r^ is relatively very small 
and since the auxiliary filament is approximated by a one- 
dimensional bar, the r-dependence of will be neglected and 
it will be represented by its value at r =0. The stress in the 
auxiliary filament may easily be obtained from (41) as 

o (z) = - p - / Z(t)dt . (50) 

^ 0 z 
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Thus, after some manipulations the axial stress 1n the filament 
is found to be 

c c 

Of (z) = On - ptij(z) - 7 ^ / z(t)dt + / h2(z,t)Z(t)dt , 

0 z -c 

(0<z<c), (51) 

where 


" 4(lU) + 


c-z 


c + z 


[r^ + (c-z)2]^/2 [r^ + (c+z)2]1/2 


(c-z) 


Lc±zl 


] . 


[r^ + (c-z)^]^^^ [r^ + (c+z)^]^"^^ 

r.(t-z) \2 

h-(z,t) = 5 — 5 -]. (52. a. b) 

2 4(l-v)[r2 + (t-z)2]3/2 r‘+(t-z)2 


The results of a numerical example giving the filament stress 
are shown in Figures 5 and 6. Figure 5 shows various 

combinations of c/r^ and E^/E. 

For large c/r^ ratios it is reasonable to expect that the 
relative contribution of the end tractions p (particularly away 
from the ends) would be negligible. Figure 6 shows the results 
obtained with and without ignoring the effect of p for various 
combinations of c/r^ and E^/E. It is clear from the figure that, 
in future calculations regarding the application of the filament 
model developed , in this paper, the effect of the end tractions 
may indeed be ignored. 

4. COMPARISON WITH OTHER MODELS 

Two other possible models for an elastic filament are the 
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Figure 5. Axial stress in an elastic filament. 




Distance From Center-z/c 


Figure 6. The effect of end tractions p on the axial 
stress in an elastic filament. 



ellipsoidal inclusion considered in [1] and the model discussed 
in [2]. The solution given in [1] is in closed form where it is 
shown that the stress state in the inclusion is uniform. The 
expression for the stresses are rather lengthy and will not be 
presented in this paper* ** . The calculated results for the stresses 
= Oqq and in the inclusion (filament) which is in the form 
of an ellipsoid with the semi-axes (cjr^jr^) are shown in Figures 
7 and 8. Figure 9 shows the comparison of the maximum filament 
stresses obtained from the ellipsoidal inclusion solution 

and from the elastic filament model given in the previous section 
(equation (51)). The agreement appears to be quite good for 
lower values of E^/E and acceptable for higher E^/E. 

Following the procedure outlined in [2] the filament problem 
considered in this paper may be reduced to a Fredholm- type inte- 
gral equation of the second kind with a logarithmic singularity. 

In this model the body force Z(z) acting on the matrix is assumed 
to be distributed uniformly over the cross-section (z = constant, 
0<r<r^) and the integral equation is obtained by matching the 
strains in the matrix and in the auxiliary filament (which 

is also assumed to be a one-dimensional bar). The calculated 
filament stresses obtained from this model for various combina- 
tions of c/r^ and E^/E are shown in Figure 10. Figure 11 shows 
the comparison of the filament stresses obtained from the models 
given in [1] and [2], and from that described in this paper. 

*The details may be found in [9]. 

**The details of the derivation of the integral equation and the 
solution may be found in [9]. 
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Figure 7. The axial stress in an elastic ellipsoidal 
inclusion [Ref. 1], 







Figure 10. Axial stress in an elastic Inclusion calculated 
from the model of Ref. 2. 
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Figure 11. Comparison of the axial stresses obtained from 
three different filament models. 


5. 


THE CASE OF MULTIPLE FILAMENTS 


The application of the filament model developed in Section 
3 of this paper to problems involving multiple filaments is 
straightforward provided the distance between any two filaments 
Is sufficiently large compared to r^ , (1 = 1 , . . . iN) (r - being the 
radius of the ith filament) so that the variation of the body 
force Z.(r^,z), (i=l,....N) along the circumference of the 
filament may be neglected. In this case the problem may easily 
be shown to reduce to a system of N singular integral equations 
in the unknown functions Z-(z), (i=l»..»»N) which may be solved 
numerically in a routine way [7*8]. In particular, if the fila- 
ments are identical and are located symmetrically the problem 
may be simplified considerably. This is the case where 

= . . . = = . . . = C(^ - c. z = 0 is a plane of symmetry for 

all filaments (which are parallel to the z axis), the filaments 
are evenly spaced on a circle of radius b on z=0 plane, and the 
matrix is again subjected to a uniaxial stress parallel to 
the z axis and away from the filament region (see the insert in 
Figure 12). For this special case because of symmetry 
Z^(z) = ... =Z^(z) =Z(z) and p^ = . . . = p,^ ^ p , and hence the 
problem reduces to a single integral equation and a single alge- 
braic equation in Z and p. Referring to Section 3 of this paper, 
after some simple manipulations we find 

J + / m2-(z,t)Z(t)dt + „ - E) ^ Z(t)dt 

ce. O' f ' Z 

2C E 

+ I [tn 2 (z) + E-Vf] = - > (-c<z<c), 
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2C,E c 


/ m,(t)Z(t)dt + - E) I *^(*>'^* 


CiE 


T 


(53. a, b) 


m 


m 


j2(z.t) = k22(z.t) + Trr^(t-z) 

X [1 . 2y + 2 ^ • 

d? + (t-z)^ 

1 

j(z| . kj(z) . .r^ ,y[<,2.,'''|!]371 <’ ■'’'*4 ’.‘u!!)2’ 

. C + Z n V + 3y(c+z)^ \-i 

[d^ + (c-z)2]^''^ d^+(c+z)^ 

,(z) . k,(z) . ,L'[d;.(c^)2]''^ " * ■ 


2 

irro N ^ y 

*^3 ‘ ^2 ■ 77 i=2V+d?/c2)1/2 - 77777' ■ '■i 


(1 + - #-] . 
i 


r 47 t(1-v) 

" (1+vj (3 - 4v[ ’ 

d^ = b{2[1 - co$(2ir{i-l}/N)]}^''^ . (i = 2 N), (54.a-f) 

vyhere ^q»c are the dimensions and is the Young's modulus of 
the filaments, E and v are the elastic constants of the matrix, 
k^(z) and l< 2 (z) are given by (47), and y are given by (44), 
l< 22 {z>t) is given by (37. c), Z(z) is the unknown body force along 
the filament-matrix interface ('''= 1 ^ 0 * |z|<c), and p is the uni- 
formly distributed body force applied to the ends (0<_r<rQ, 

2 = + c) . 
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Some of the numerical results obtained from the solution of 
(53) are given by Figures 12-14. Figure 12 shows the distribu- 
tion of the axial stress in two symmetrically located filaments, 
where again the problem is solved with and without taking the end 
effect p into account. In this example, the effect of p is seen 
to be quite insignificant. Again, for N=2, the effect of length- 
to-diameter ratio and spacing of the filaments are shown in 
Figure 13. Note that as the distance b between the filaments 
decreases the axial stress in the filaments also decreases. 

Figure 14 shows the ratio of the maximum filament stress (which 
is at z=0) for N>1 to that for N=1 as a function of the distance 
parameter b (see insert in Figure 12) (Figure 14. a), and as a 
function of modulus ratio E^/E (Figure 14. b). As expected, the 
interaction effect increases with increasing N and increasing 


E^/E, 
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Figure 12, The axial stress in two identical filaments obtained 
with and without including the effect of the end tractions p. 
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Fiaure 14 The ratio of the maximum filament stress for N>1 
ti'that for N= 1 as a function of the distance parameter b 
and the modulus ratio E^/E. 
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APPENDIX A 


Integrals used in the derivation of the kernels 
(i , j - 1 .2) given by (32) : 


= 2r^(l - cosa) + (t-z)^ , 

T— =• , , 1/j . = 

0 P [4r2+ (t-z)2]''" 


4r: 


4r^ + 


cosada _ 4K(k) + ( t-z ) _ ) 

. 2 [Ir* . 

^0 

Y ^ - 4E(kJ ^ , 

0 (t-z)^[4r2 + (t-z)^]'^" 

(1 - cosa)^ , 2 [K{k) - E{k)] 

0 r^[4r2 + (t-z)2]l/2 


n-cosa)^da . 2[2r^(t-z)^] 

I p2 r^[4r2 + (t-z)2]^/2 


E(k) 


r^[4r2 + (t-z)^]’^'^ 
0 0 


(A.l) 
(A. 2) 


(A. 3) 


(A. 4) 


(A. 5) 


(A. 6) 
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APPENDIX B 


Simplified expressions for the functions M^.(z) and N^.(z) 
given by (45): 

Note that 

dN. (z) 

— ^ - (2i -l)(c-z)N^^^(z) . 
dM. (z) 

— — = - (21 - 1 ) (c+z)M^^^ (z) . 


(B . 1 a , b) 


A1 so 


2tt 

/ 


de 


4K(ki) 

2 


0 [r^ + r^ - 2rr^cos0 + (c±z)^] [(r+r. + (c±z)^] 
7 4rr_ 


27T72 


, (B.2) 


1 ' (Vr)2Mc±z)2 ’ 

where the upper and lower indexes in the modulus k correspond to 
the upper and lower signs in (c±z). Thus, 

r« 


rK(k,)dr 

M-i(z) “ 4 J g i“y p » 

^ 0 [(r+r^)^ + (c+z)2] '' 

•"o rK(k-)dr 

Ni(z) = 4 / s — ^ TTTl > 

^ 0 [(r+r„)^ + (c-z)^] ' 


(B.3a,b) 


and 


M2(z) 
N2(z), 
M3(z) 


= 4 / 


rE(ki )dr 
2 


8 [{r+r^,)2 + (c±z)2]l/2[-(r^.^)Z +{c±z)2] 


L a 4 I 

Nq(z) J ^0 [(r +r)^ + (c±z)^]^^^[(rQ-r)^ + (c±z)^] 


rdr 


-K(ki) 

2 


(rQ+r)*^ + (c±z) 

1 

(r^-r)^ + (c±z)^ 


^ + 2E (k] ) [ 


2 (I'n+r)^ + (c±z)^ 


]} ■ 


(B.4a,b) 
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